Unitarity violation in sequential neutrino mixing in a model of 

extra dimensions 



Subhaditya Bhattacharya 1 , Paramita Dey 2 and Biswarup Mukhopadhyaya 3 
Regional Centre for Accelerator-based Particle Physics 
Harish- Chandra Research Institute 
Chhatnag Road, Jhunsi, Allahabad - 211 019, India 



Abstract 

We investigate the possibility of unitarity violation in the sequential neutrino mixing 
matrix in a scenario with extra compact spacelike dimensions. Gauge singlet neutrinos 
are assumed to propagate in one extra dimension, giving rise to an infinite tower 
of states in the effective four-dimensional theory. It is shown that this leads to small 
lepton-number violating entries in the neutrino mass matrix, which can violate unitarity 
on the order of one per cent. 
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1 Introduction 



The ever-consolidating evidence in favour of neutrino masses and mixing has spawned a 
large volume of speculations on new physics possibilities that could be at their origin. Con- 
sidering the three light sequential neutrinos, many proposed scenarios, including seesaw 
models of type I [1], II [2] or III [3], ensure unitarity to a high degree of precision in the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix describing mixing in the lepton sector 
[4]. A measured departure from such unitarity, evinced from precision data in the neutrino 
sector, may thus point towards some novel mechanism for the generation of neutrino masses. 
One such possibility arises when at least one small gauge singlet Majorana mass term enters 
into an extended neutrino mass matrix. Under certain conditions, this situation passes off as 
'inverse seesaw mechanism' [5, 6, 7]. It has been demonstrated in a number of recent works 
that this can lead to a violation of unitarity at the level of about 1% or more in the 3x3 
light sequential neutrino mass matrix, due to mixing with additional sterile states [5, 6, 7, 8]. 
A pertinent question to ask is: does such a situation fit into some of the popular scenarios 
of new physics at the TeV scale? 

The experimental constraints on the loss of unitarity as well as its testability in neutrino 
oscillation experiments has been investigated recently [8, 9]. As for theoretical models, a 
GUT-inspired scenario, based on SO(10) with a breaking chain involving an extra U(l) 
gauge symmetry surviving at low scale, has been considered recently for this purpose [6]. 
This scenario has been shown to lead to two-loop generation of some small Majorana masses 
and consequently lead to the inverse seesaw mechanism. It has also been suggested that a 
super symmetric model including two types of gauge singlet neutrino superfields may produce 
effects of this kind [10]. In a number of other model-building ventures, too, the effect 
mentioned above emerges as a consequence [11]. 

Phenomenological implications of unitarity violation in the PMNS matrix, including its 
signatures in phenomena driven by neutrino oscillation, have been recently investigated [9]. 
In this paper, we point out that a loss of unitarity in the PMNS matrix can also arise if one 
has extra flat spacelike dimensions, with gauge singlet neutrinos propagating in one extra 
dimension, and lending small diagonal elements to an extended neutrino mass matrix. 

Mechanisms of neutrino mass generation have been frequently suggested in models of 
compact extra spacelike dimensions, both flat [12, 13, 14, 15, 16, 17, 18] and warped [19]. 
Here we consider a minimal higher-dimensional framework where the standard model (SM) 
fields all lie on a 3-brane, while one or more gauge singlet neutrino propagate along one 
flat extra dimension [12, 13]. However, there can in principle be several extra spacelike 
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dimensions where gravity propagates, thereby evading the already established lower limits 
on the number of such dimensions [20]. An orbifold symmetry is further imposed along the 
compact direction containing the neutrino(s), so that one obtains only one (right-handed) 
chirality for the n = Kaluza-Klein (KK) mode. It has been demonstrated earlier that this 
scenario can naturally suppress neutrino masses via a Type I [1] seesaw mechanism. 

The gauge singlet neutrinos can have Majorana masses in five dimensions to start with. 
We are specially interested in the situation where a strong cancellation between this mass and 
the KK tower mass leads to a very small entry in the effective neutrino mass matrix in four 
dimensions. We show that the resulting mass matrix has additional 'sterile' states mixing 
appreciably with the sequential neutrinos. It is found that one can consequently expect the 
violation of unitarity in the 3x3 (PMNS) matrix in certain regions of the parameter space 
of such a model. 

In Section 2, we outline some scenarios that lead to loss of unitarity of the PMNS matrix, 
including the inverse seesaw mechanism. The extra-dimensional model under investigation 
is briefly reviewed in Section 3. The viability of a substantial loss of PMNS unitarity is 
numerically demonstrated in Section 4. We summarise and conclude in Section 5. 

2 Loss of unitarity in the PMNS matrix 

In general, the well-known Type I seesaw mechanism involving a light and a heavy neutrino 
also involves a departure from unitarity in the PMNS matrix. However, this departure is 
immeasurably tiny, since the seesaw mass scale in invariably much higher than the light 
neutrino masses [1]. An exception to this may occur if a very small Majorana mass is 
introduced. However, this choice is inhibited by (a) the need of justifying such a small 
AL = 2 mass in terms of new physics, and (b) the need of introducing excessively suppressed 
Dirac masses for generating light neutrinos, which essentially destroys the motivation of the 
seesaw mechanism. 

The situation can be different when one has more than one two-component sterile neutri- 
nos. It has been shown in a number of recent works [5, 6, 7, 8] that this allows one to insert 
small lepton-number violating mass terms in diagonal positions of the neutrino mass matrix, 
and the off-diagonal entries need not all be much smaller. Its most noticeable consequence 
is a loss of unitarity at the level of 1% or more in the 3x3 (PMNS) part of the neutrino 
mass matrix. 

Several kinds of scenarios that meet this description are found in the literature [5, 7, 8]. 
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Here we mention two classes only among these, considering for illustration one sequential 
and two sterile species in each case. The first [5] is one of the form 



M = 



1 m D \ 

m D m R m N 
\ m N m L J 



(1) 



in the basis {pi,, N r , Nl), where the last two are gauge singlets. The masses rriL t R arise from 
AL = 2 terms. For << mp, m^, this not only yields an active neutrino mass eigenstate 
in the right order, but also leads to active-sterile mixing at the level of 1% for appropriate 
choice of the mass parameters (say, for example, mo ~ 10 GeV, ~ 1 TeV, and mi ~ 10 
keV) . A corresponding situation with three sequential neutrinos will show unitarity violation 
at the same level in the PMNS matrix, but with an additional light sterile neutrino. Since the 
light (sequential) neutrino mass vanishes in the limit m L — > 0, it is often called an 'inverse 
seesaw' scenario. 

Another situation that one can consider has the same choice of neutrino basis states, but 
a mass matrix of the form [8] 



M 



\ 



' m D m N 

m D m R 
\m N m L ) 



(2) 



with mo « rriR and m 2 N /rriR « « « rriR. It has been found that this 
situation, too, leads to light sequential neutrino(s) and unitarity violation at the same level 
(~ 1 %), for appropriate choice of parameters (say, for example, m D ,m N ~ 1 MeV, m R ~ 1 
TeV and m L ~ 100 eV). The difference with the previous situation is that (a) one obtains a 
light sterile neutrino even with one sequential family, and (b) the sequential neutrino mass 
does not vanish in the limit mi — > 0. This makes it deviate from an inverse seesaw scenario 
in the strict sense, although it is equally interesting from the viewpoint of unitarity loss 
of the PMNS matrix. Since such unitarity loss is a very interesting consequence that is 
experimentally testable, it is worth exploring if it occurs in some otherwise well-motivated 
theories beyond the SM. In the next two sections we outline one such scenario, and go on to 
examine its potential for generating unitarity loss. 
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3 A model with extra dimensions 



In this section we describe the model adopted for illustrating our point. It assumes extra flat 
compact spacelike dimensions where gravity can propagate. The SM fields are confined to a 
3-brane which constitutes a 'slice' in the higher- dimensional space. So far it is very similar 
to the Arkani Hamed-Dimopoulos-Dvali (ADD) scenario [21], excepting that it includes an 
effort to account for neutrino masses, through the introduction of one gauge singlet neutrino 
propagating in one extra dimension only [12, 13, 14, 15, 16, 17]. Thus, while all the phe- 
nomenology involving gravitons remains similar to that in the ADD framework with several 
extra dimensions, one can consider just the five-dimensional subspace for studying neutrino 
physics. We make our analysis simple by adhering to one generation of SM neutrinos. The 
fifth flat dimension (y), along which propagates the right handed neutrino (N(x,y)), is com- 
pactified over an S 1 /Z 2 orbifold where R is the radius of compactification. The preservation 
of the Z 2 invariance necessitates the existence of at least two symmetrically placed branes, 
and the SM fields lie on either of them. Thus the complete leptonic field content of the 
model is 

\£l(x) J \r]{x,y) J 

where V£, £l, £r are Weyl spinors in four dimensions, and £, rj are two-component spinors in 
five dimensions. Under S 1 /Z 2 , the latter may be associated with opposite parities: 

t(x,y) = £(x,-y), v(x,y) = -rj(x,-y). (4) 

The brane where the SM is localised, can be assumed to be at y = a just for generality, 
instead of at the orbifold fixed point y — 0. We shall see later that this adds to the freedom 
of the model. The generic effective four-dimensional Lagrangian of this model is given by, 

2vr_R 

C eS = Jdy {iV(V«9 M + l5 d y )N - - (MN T C^~ l l5 N + h.c.) 





+ S(y-a) 



+ T^yf 2 L&V + h.c. ] + 5(y- a) C SM } , (5) 



L(M F )V2 s {Mp) 



where $ = ia 2 Q*, Csm is the SM Lagrangian, M is the Majorana mass for iV (we do not 
specify its scale for the moment), is the 5-dimensional charge conjugation operator and 
Mp is the fundamental gravity scale. The Yukawa couplings in five dimensions, h± t 2, are 
assumed to be 0(1). For gravity propagating in a d- dimensional bulk, 

M P = (2irM F R) d/2 M F , (6) 
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for the simple case where all the compactification radii are of equal size R, Mp being the 
four-dimensional Planck scale. A Dirac mass term m D NN is not allowed in equation (5) 
because of the Z 2 symmetry. 

Following equation (4), the two-component spinors £ and 77 can be expanded as, 



rj(x,y) 



1 



V2tiR 
1 



1 



-kR n=1 
ny 



E ^(x) COS ( ^ 



'txR 



J2 Vn(x) sin ( — 



n=l 



(7) 
(8) 



where the chiral spinors £ n (x) and r) n (x) form an infinite tower of KK fields. Using these 
expansions and integrating out the y-coordinate, the effective Lagrangian reduces to 



£ cff = C SM + £0(^)6 + (h?L$*Z - -M£ £o + h.c.) + £ 

1 71=1 

+ fj n (ia IX d ll )r) n + ^(^nVn + £nVn) ~ ^ M (f n f n + + h.C.) 

+ v / 2(M n) L$*£ n + h^L^rjn + h.c.) 
in a basis in which M is positive, and with: 



hi 



(n) 



(n) 



(2tt MpR) 1 ! 2 
h 2 

(2irM F Ry/i 



cos 



sin 



na 
~R 
na 
~R 



M F \ 1 l d 



M P 
M F 



h\ cos — 



l/d 



h 2 sin 



na 
R 
na 
~R 



hi cos 
h 2 sin 



£„(^^)£„ 

(9) 

(10) 
(11) 



na 
na 



For deriving the last two equalities on the right hand sides of equations (10)-(11), we have 
made use of equation (6). 

Equations (10) and equation (11) imply that the induced four-dimensional Yukawa cou- 
plings J112 can § e t suppressed by many orders depending on the hierarchy between Mp and 
Mp; for example, if gravity and the bulk neutrino feel the same number of extra dimensions, 
say d — 1, then these couplings are suppressed by a factor M F /M P ~ 10~ 15 , for M F m 10 TeV 
(see also [12, 13]). 

It is clear from equation (3) that £ and fj have the same lepton number. Thus, the 
simultaneous presence of the two operators and L§*n in equation (9) leads to lepton 

number violation. Such coexistence of the two operators is possible only if we allow the brane 
to be shifted by an amount 0) from the orbifold fixed points (y = 0, ttR). Such a shifting 
of the brane, respecting the Z 2 invariance of the original higher dimensional Lagrangian, 
has been shown to be possible under certain restrictions in Type-I string theories [22]. As 
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indicated in [12, 18], the Z 2 invariance can be taken care of by allowing the replacements 

1 



r)8(y - a) 



2 
1 



£ 8(y -a) + 5(y + a- 2nR) 



rj 



5(y - a) - 5(y + a- 2nR) 



(12) 



with < a < 7[R and < y < 2nR. Here we re-iterate that a Z 2 -invariant implementation 
of brane-shifted couplings requires the existence of at least two branes placed at y = a and 
y = 2nR — a. 

A remarkable feature of the brane-shifted framework was pointed out in [18], where it 
has been shown that in such a framework it is possible to completely decouple the effective 
Majorana- neutrino mass (to) and the scale of light neutrino masses, so as to have (to) within 
an observable range. Therefore, the Lagrangian (9) contains two types of Majorana neutrino 
mass terms (involving respectively the parameters M and h^) both of which lead to a 
breaking of L. Such L-breaking is a necessary ingredient of leptogenesis. 

Following the notations of reference [12], we now introduce the weak basis for the KK 
Weyl-spinors, by defining 

X±n = ^= (£n ± Vn), (13) 

followed by a rearrangement of the states £ and xt.-* sucn that, for a given value of n (say, 
n = k ), the smallest diagonal entry of the neutrino mass matrix is 



mm 



( 



M — — 
R 



) < l/(2i2). (14) 
After re-ordering, we can define the multiplet ^ u consisting of the Majorana spinors 









while the effective Lagrangian for right handed neutrinos reduces to 



(15) 



(16) 
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where A4^ K is the corresponding neutrino mass matrix given by 



M* K = 



( 


m (0) 




1) 


m (!) 




2) 


m (2) 


m (0) 


£ 


















mi- 1 ) 





£ — 


1 

R 













m (!) 










Z+R- 























e — 


2 
R 





m (2) 





















V 



(17) 



The most important consequence of such a rearrangement is that the mass scale M, which 
we did not specify earlier but which could be arbitrarily large, is now replaced by the light 
mass scale e. The entries in the first row and the first column of M^ K are given by the 
relation, 

(n — A; )a N 



m 



with 



(n) 



V 

72 



hi cos 



(n — k )a 
R 



+ h 2 sin 



R 



m cos 



na 
~R 



(18) 



v hf + h 2 2 
m — — 



2 V 7iM F R y/M^R' 
4> h = tan" [ Vi )+k -, 

where v is the vacuum expectation value of the SM Higgs boson. 



(19) 
(20) 



4 Unitarity loss with extra dimensions: some numeri- 
cal illustrations 

Here we show that a substantial loss of unitarity of the PMNS matrix can occur in different 
allowed regions of the parameter space of the model described in the previous section. The 
first issue is, of course, insuring at least one small entry in diagonal positions of the neutrino 
mass matrix M^ K . Equation (14) tells us that e < 1/2R. There is no other theoretical or 
phenomenological constraint on e. Thus e qualifies to be the small diagonal element which 
can be potentially responsible for a departure from unitarity. 

Mp, the Planck mass in five dimensions, is expected to be ^ TeV, since gravitational 
effects will otherwise become important in low-energy physics. At the same time, in order 
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to ensure that physics along the compact dimension(s) is not plagued with trans-Planckian 
effects, one should have l/R < M F . Thus, in the expression for m in equation (19), M F R 
is at least of order unity. Given the fact that the five-dimensional Yukawa couplings hi t2 , 
too, are prima facie of the order of unity, this implies that m can at most be around v, the 
electroweak symmetry breaking scale. 

We indicated two scenarios of unitarity loss in section 2: (i) where the off-diagonal 
elements in the extended neutrino mass matrix are all larger than the diagonal ones, and 
(ii) where diagonal elements excepting the smallest one are larger than the off-diagonal ones. 
Let us first examine whether the extra-dimensional model under scrutiny answers to both of 
these scenarios. 

The first possibility demands 

m > \ (21) 

since, with a ^ 0, can approach m for some value of n along the tower. Using equation 
(19) in (21), one obtains 

max ^ _ fOO\ 

Mp R [ZZ) 
The inequality should hold for the maximum value of the right-hand side for a given Mp, 

which is M F itself. Thus we have 

M F < m max (23) 

Therefore, demanding m > l/R implies that the five- dimensional Planck scale has to be 
brought down below m max which is just about the electroweak symmetry breaking scale, and 
hence is inadmissible. 

The first scenario is thus disfavoured in this model. On the other hand, since the diagonal 
elements M^ K (i, i), % > 3 are always greater in magnitude than for all values of n, one 
can say that for sufficiently small e, this model provides an opportunity for unitarity loss in 
the PMNS matrix in the sense of the second scenario mentioned in Section 2. 

The value of l/R, on the other hand, is not subject to any general constraint stronger than 
that arising from the validity of Newton's law of gravitation down to about lCT 2 mm, which 
essentially allows l/R to have any value ^ 1CT 8 MeV [20]. Precision electroweak constraints 
do not tighten the constraint, since the tower resulting from the compactification of the 
extra dimension corresponds to gauge singlet neutrinos only. Thus, in order to have loss of 
unitarity in the PMNS matrix, we are faced with two possibilities, namely, (a) e « l/R, 
and (b) e ~ l/R. We show below that both of these situations are possible. 

Our principal aim is to check if it is possible to have the sequential neutrino masses in 
the right order (~ 10~ 2 eV), and at the same time have substantial violation of unitarity. 
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The latter requires that the squares of elements of some particular column in the full mixing 
matrix beyond the PMNS block add up to (9(1CT 4 ). This sum is defined as 5 2 here. 1 

We carry out this investigation in the simplified situation, with one sequential neutrino 
flavour and just one gauge singlet neutrino in the bulk. We shall comment later on the 
generalisations necessary to generate the actual pattern of masses and mixing. As far as the 
violation of PMNS unitarity is concerned, however, the conclusions we reach below remain 
valid even when such generalisations are made. 

Case (a): This implies a fine cancellation between the bulk mass M and some integral 
multiple oil/R. While a dynamical explanation of this is difficult, it is not entirely unlikely, 
as both M and 1/R can rather naturally be around the TeV-scale, and there is a distinct 
possibility of the two of them having near-coincident values. 

A few illustrative points in the parameter space for this case are shown in Table 1. We 
have confined ourselves to 1 MeV < 1/R < 10 TeV. It is found that, to get a substantial 
unitarity violation (5 > 0.5 %) and neutrino mass in the right order, the largest possible 
value of e that we can take is ~ 10~ 6 . In this case, the 2x2 block in the upper left corner of 
the neutrino mass matrix effectively determines the masses of the sequential and the lightest 
sterile neutrino, given as (W°)) 2 / 'e and e respectively. One further has: 

^1} — (24) 

Therefore if e is increased, the concomitant enhancement in required to keep the 

sequential neutrino mass unaffected, ends up suppressing 5. The values of required 
point towards 4>h ~ 7r/2. On the other hand, the fact that rrS 1 ^ can vary over a wide range, 
implies that the brane-shift parameter a can vary from zero to 0.1R approximately. It should 
be noted that Table 1 includes one sample corresponding to = mP-\ which means a =0. 
Thus, in this case, large unitarity violation is consistent with both the cases where the brane 
is located at orbifold fixed point and where it is noticeably shifted. 

Case (b): In this case e and 1/R can be relatively close to each other. Therefore, no 
drastic cancellation is required between them, and no allegation of fine-tuning can be levelled 
against such a scenario. 

Some sample results for this case are presented in Table 2. It should be noted that 
values of 1/R of comparable smallness as that of e implies that the part of Mf K beyond 
the upper left 2x2 block no longer tends to decouple. An immediate consequence is that 



1 In references [8, 9], the violation of unitarity has been defined in terms of a parameter r\. It is easy to 
check that 2?7 = S 2 . 
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Table 1: Different sample points in the parameter space of the model where substantial unitarity 
violation takes place, for e « 1/R. All mass parameters are in GeV. 

m^ 1 ^ have to be as small as when 1/R is small. This in turn drives the value of the 
brane-shift parameter a to very small values. Therefore, unlike in the previous case, the 
brane is compelled to be close to the orbifold fixed points. 

The figures included in Table 2 are self-explanatory; the provision for substantial unitarity 
violation is clearly there. However, there is a relative paucity of available points compared to 
the previous case. This is because the part of Ai^ K beyond the 2x2 block is not ineffectual 
in determining the sequential neutrino mass and its mixing with light sterile states. In order 
to comply with all constraints there, one therefore requires a correlation between e and 1/R, 
in contrast to the situation with e << 1/R, thereby restricting the allowed points in the 
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parameter space. 

The numerical results presented by us are obtained through the diagonalisation of the 
4x4 neutrino mass matrix, including a tower of states up to the first KK excitation only. 
We have, however, checked that the results do not change qualitatively upon the inclusion 
of additional towers and the resulting augmentation of the mass matrix. 



l/R 
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m (0) 




H%) 
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10" 8 


3.5 x IO" 8 
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KT 5 


2 x 1(T 6 


1(T 8 
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5 x 1(T 7 
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Table 2: Different sample points in the parameter space of the model where substantial unitarity 
violation takes place, for e c± l/R. All mass parameters are in GeV. 



5 Summary and conclusions 

We have studied a popular model of flat extra compact spacelike dimensions. A gauge singlet 
neutrino is assumed to propagate in one extra dimension. The 'bulk' mass possessed by this 
neutrino can undergo cancellation with the KK tower mass for some member of the tower, 
giving rise to at least one small diagonal entry (e) in the infinite-dimensional neutrino mass 
matrix in four dimensions. We show that this can cause substantial mixing between the 
sequential and sterile neutrinos without violating any existing constraint. The consequence 
is a departure from unitarity of the PMNS matrix, both for e « l/R and e ~ l/R, R being 
the radius of the compact dimension housing the gauge singlet neutrino. 

For the sake of simplicity, we have presented our results for one sequential neutrino. It 
can be easily checked that the conclusions are valid with additional generations. In fact, the 
constraints on 5 are easy to satisfy, since the strongest constraint on the PMNS matrix is 
on its (1, 2)Hh element [8, 9]. The relatively unconstrained mixing of, for example, v T with a 
sterile neutrino can accommodate the values of 5 obtained here. On the other hand, it may 
be difficult to accommodate the neutrino mixing data and mass hierarchies with one sterile 
bulk neutrino only [23]. At least two such neutrinos can, however, accommodate everything 
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rather easily, thanks to the additional Yukawa couplings available, which are essentially free 
parameters. Our general conclusions are unaffected by such extensions. 

In conclusion, the phenomenon of unitarity violation in the PMNS matrix can be moti- 
vated rather well in a model of extra dimensions. This brings to the fore the likely connection 
between subtleties of the neutrino sector and theories which advocate strikingly new physics 
around the Tev scale. 
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